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Abstract. In the recent years we have learned that light quarks play a 
crucial role in QCD-like theories, transforming it to many different phases. We 
review what is known about them, both from lattice and non-lattice approaches. 
A particularly simple mechanism of the QCD chiral restoration phase transition is 
discussed first: it suggests that it is a transition from randomly placed tunneling 
events (instantons) at low T to strongly localized tunneling-anti-tunneling pairs 
at high T. Many features of the transition found on the lattice can be explained 
in this simple picture. Very relevant for RHIC, this approach predicts a strong 
non-perturbative interaction between quarks above the phase transition. It also 
predicts that QGP-like phase sets in at zero temperature, provided few more light 
quark flavors are added to QCD. Finally, we also discuss possible experimental 
signatures of the QCD phase transition. One issue is CERN dilepton data, possibly 
related with "dropping" masses of p, A\ mesons. Another is direct manifestation of 
a softeness of EOS (smallness of pressure/energy density) in the phase transition 
region in flow and even the global lifetime of the system. 



1 Introduction 

My topic is to review our current understanding of critical phenomena in 
QCD (and its close relatives). But before that, let me comment on a some- 
what different type of critical phenomena discussed in one (late) session of 
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the workshop, the so called "self-organized criticality" exemplified by the fa- 
mous sand piles. Not being an expert in this field, I still dare to suggest that 
the RHIC summer workshop itself is a perfect example of this phenomenon. 
Indeed, only several years ago many of the participant of RHIC project (ex- 
perimentalists, accelerator people and even most of the theorists) would not 
even listen to a talk with such title as mine. (This is commonly described 
by exponentially decaying correlations, with rather short correlation length.) 
Now, with common goals and concerns about directions future experiments 
at RHIC will follow, we are all well inside "one correlation length". This is 
clearly a phase transition in its own right. 

Returning to QCD, let me start with the comment that during the last 
few years we have learned about many new phases, which the gauge theories 
may have. Of course, the one we are going to look for experimentally is still 
the chirally symmetric Quark-Gluon Plasma (QGP), in which the charge is 
screened |]J rather than confined. Recent unexpected findings suggest that 
it may have completely new features, such as preserve some hadronic modes 
as a bound states . 

However, theory predicts some other phases which has taught us many 
new lessons. In particular, it seems that QGP may exist even at zero temper- 
ature in QCD with about 5-7 light fermions. Even more flavored QCD, with 
7-16 light quarks, is expected to be in even stranger conformal phase, an 
extended relative of a condition commonly studied only at the second order 
phase transition points. Similar phases are proven to exist in supersymmetric 
extension of QCD, in which there were recently a significant advances due to 
Seiberg and Witten. 

In addition, we have learned about some unwanted phases, such as (i) the 
Aoki phase appearing in lattice studies with Wilson fermions (see Ukawa re- 
view) and (ii) the Stephanov phase appearing in quenched lattice simulations 
with the non-zero chemical potential. 

As emphasized by T.D.Lee in his opening talk, before thinking about 
"small" (the particles) one should first clearly understand "large" (the phase 
one actually is in): as we will see it is a very good advice to these lattice 
calculations indeed. 

The bottom line of this talk is the crucial (but still not quite well under- 
stood) role of light fermions, generating all these phases of QCD. In particu- 
lar, compare the phase transition found in "quenched" calculations (Nf = 0) 
and those with Nf = 2 — 3 dynamical quarks. The former case has a "de- 
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confinement" phase transition at high T c « 260MeV , while the latter show 
a "chiral restoration" transition at much lower T c « ISOMeU [] Unlike at 
normal(T=0) conditions, at T « T c physical quantities are very sensitive to 
such little details as the mass value of the strange quark. Therefore we are 
still not quite sure about the order of the transition: the latest lattice results 
[0] incline again toward the 1-st order in the real world. 

But do we understand why the transitions happen? Is there a simple 
picture which can explain its microscopic mechanism? Can we build a work- 
ing model, reliable enough to provide some guidance in delicate questions 
relevant to experimental observables? 

Qualitative explanations of why the QCD phase transition takes place are 
often done in an over-simplified way, emphasizing the "overlapping" hadrons 
in a bag-model-type picture. But such pictures give all numbers f\ and physics 
wrong. Pure gluodynamics is an especially good example: at T c fa 260MeV 
the density of glueballs is negligible since even the lightest one has a mass of 
about 1.7 GeV! Looking at the low-T side of this transition it is impossible to 
tell why it happens. The same is qualitatively true for the chiral restoration: 
on the hadronic side the matter is still relatively dilute. But the reason why 
the phase transition happens is very simple: just two very different phases 
happen to have the same free energy. The lesson: one cannot understand a 
QCD phase transition without a quantitative model for both phases. 

The energetics of the "deconfinement" and "chiral restoration" transitions 
is entirely different, suggesting different physics. The former has huge latent 
heat, few GeV/ fm 3 , so that probably all of the non-perturbative vacuum 
energy density (proportional to the "gluon condensate") is "melted" in it. 
This in not the case for "chiral restoration" in QCD: large portion of the 
gluon condensate should actually survive it U. What is this remaining "hard 
glue" or "epoxy", as Gerry Brown called it? Why, unlike "soft glue" does 
it not produce a quark condensate? Can it affect quark interactions and 
hadronic masses? 

The major "non-lattice" approaches to the problem include: (i) mod- 
els based on particular effective Lagrangians, like the sigma model or chiral 

1 By changing the quark masses from light to heavy continuously, it was shown that 
these two transitions are indeed different phenomena, separated by a large gap in which 
there is no transition at all. 

2 For example, the MIT bag model literally predicts QGP formation in heavy ion colli- 
sions at unrealistically small (BEVALAC/GSI) energies. 
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effective Lagrangians models |5|, [J; (ii) QCD sum rules at finite tempera- 
ture/density 0; (iii) the interacting instanton liquid model (IILM) |2], ||. 

Effective Lagrangians are very powerful tools at low temperatures: with 
parameters fixed by data one can indeed accurately account for effects due to 
non-zero occupation factors. However, they are clearly "one-sided" , unable to 
deal with QGP. In principle, QCD sum rules are not limited to a description 
of only hadronic phase: the structural changes can be adequately described 
by VEV of different operators, or "condensates" . However, since in practice 
those are unknown, people use simplifications. The most popular one is 
"vacuum dominance", reducing average values of any quark operators to 
powers of < qq >, therefore missing non-perturbative effects above T c . 

Among theory issues discussed intensely during the last couple of years 
[§, 1(| [IIJ |12j is the fate of the U(l) chiral symmetry. In cannot be exactly 



restored, as suggested in JITl, but its violation certainly is dramatically re- 
duced at T 7n T c . In practice it means that rj non _ strange (a combination of 
1], rf) and isovector scalar (we call 5) may be nearly as light as a pion. 

One very important issue (which should have been discussed more) is 
what happens with the QCD phase transition for larger number of quark 
flavors Nf. Several phase transition lines are expected there, the lowest 
probably being the chiral symmetry restoration at T=0. 



2 A mechanism for the chiral phase transi- 
tion 

The main point of this talk is that due to recent developments it becomes 
increasingly clear what is the microscopic mechanism underlying this phase 
transition. It is rearrangement of instantons, from relatively random liq- 
uid at low T to a gas of instanton anti-instanton "molecules" Q These II 
molecules are the "epoxy" mentioned above: at T > T c they indeed do not 
create the quark condensate but contribute to the gluon one. Furthermore, 
they generate new type of the inter-quark interaction, and may even create 
hadronic states, even above T c ! 

Several events lead to these developments: (i) Demonstration that "in- 

3 Note a similarity to Kosterlitz-Thouless transition in 0(2) spin model in 2 dimensions: 
there are paired topological objects (vortices) in one phase and random liquid in another. 
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stanton vacuum" explains the QCD correlation functions and hadronic spec- 



troscopy at T=0 1 13, Eflj (ii) Lattice confirmation of instanton liquid pa- 



rameters and of the dominant role in general [fH|, T5\ ; (iii) As shown in [|TB 
instantons cannot be screened^ at T < T c ; (iv) Confirmation of this statement 



on the lattice fll7| ; (v) Discovery of polarized II molecules at T rs T c [|K 
(vi) Numerical simulations and analytic studies || of the phase transition 
in the interacting instanton ensemble. 

It is very easy to explain what happens at T « T c in this approach. Re- 
call that finite temperature is described in Euclidean space-time by periodic 
boundary conditions, with the Matsubara period 1/T. So a rising T means 
a decreasing box, and when it fits to the size of one // molecule, one gets a 
"geometric" transition. In Fig. |T] from0 one can see it clearly. (The plots 
show projections of a four dimensional box into the z axis- imaginary time 
plane. Instantons and anti-instanton positions are indicated by + and — 
symbols. The lines correspond to strongest fermionic "bonds".) Notice that 
molecules are strongly "polarized" in the time direction and that they are 
separated by half Matsubara box in time At = l/(2T).f] 

In a series of recent numerical simulations it was found that like QCD, 
the instanton model has second order transition for Nf = 2 massless fla- 
vors, but a weak first order one for QCD with physical masses. Further- 
more, the thermodynamic parameters, the spectra of the Dirac operator, 
the T-dependence of the quark condensate and various susceptibilities, the 
screening masses are all consistent with available lattice data. 



3 Phases of QCD with more quark flavors 

In this section we discuss further the role of quarks in QCD, adding more 
flavors to it. If we add too much of them, namely Nf > 33/2 (here and below 
we imply 3 colors), the asymptotic freedom is lost and we get uninteresting 
field theory with a charge growing at small distances, basically a theory as 

4 Previously considered scenario based on the "instanton suppression" does not in fact 
work until 1 rather high T, well in the QGP domain, and thus it cannot be the reason for 
the phase transition. 

5 Thus one can even find the approximate phase transition temperature. The molecule 
fits onto the torus if Ap ~ 1/T C , and with known instanton size p ~ 0.35 fm one gets 
T ~ 150 MeV, close enough to the observed one. 
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Figure 1: Typical instanton configurations for temperatures T= 76 and 158 
MeV. 

bad as QED! So our "most flavored" QCD (with 16 flavors, and work out 
way down from it (see fig |H). The phase we are in there is actually rather 
simple one, known as the Banks-Zaks conformal domain P5] . It has the 
infrared fixed point at small coupling gl/l6ir 2 = —b/b' << 1 (b, b' are the one 
and two- loops coefficients of beta function). It happens in the perturbative 
domain, so the charge is small both at small and large distances. There are no 
particles in this phase, and all correlators decay as powers of the distance. In 
this phase the non-perturbative phenomena like instantons are exponentially 
suppressed, exp(— const / gl) . However, as one decrease the fermion number, 
the fixed point gl moves to larger values and eventually disappears. Lattice 
simulations of multi-flavor QCD were recently reported in |2(J. These authors 
studied QCD with up to 240 flavors. Studying the sign of the beta function 
in the weak and strong coupling domains, they confirmed the existence of an 
infrared fixed point as low as at Nf = 7. 

It is not known which phase we find next, at Nf = 5 — 6, most probably 
it is the so called Coulomb phase (basically QGP). 

The results of the interacting instanton model are summarized by Fig.|2](b- 
d) show how one singular point at Nf = 2 develops into the discontinuity line 
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Figure 2: Phase diagram (temperature-quark mass plane) of the instanton 
liquid for different numbers of quark flavors, iVj=2,3 and 5. The open squares 
indicate non-zero chiral condensate, while solid one indicate that it is zero. 
The dashed lines show the approximate location of the discontinuity line. 

for Nf = 3. The value of T c goes down with increasing Nf and at Nf = 5[] 
one finds that the chiral symmetry is restored even at T=0, provided quarks 
are light enough. 

New lattice results have been presented at this meeting by R.Mawhinney 
for QCD with Nf = 4. Details can be found in his (and N.Christ's) talk, and 
I should only say that when they have extrapolated the measured masses to 
quark masses m — > 0, they have found a dramatic significant drop in chiral 
symmetry breaking effects, such as ir — a, p—a\, N — N*(l/2~) splittings, very 
much unlike the Nf = — 3 studied before. It suggests that chiral restoration 
is nearby, very similarly to what was found in the instanton calculations. 

Finally, we return to Figure |3], and explain its rhs, showing similar phase 
diagram for the N=l SUSY QCD based on f27j. Not going into details, let 
me only mention that the vacuum is now definitely dominated by instanton- 

6 Note that Nf = 4 case is missing. It is because I have found the condensate to be small 
and comparable to finite-size effects. In order to separate those one should do calculations 
in different boxes, which is time consuming. 
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antiinstanton molecules, and their contribution can be calculated without 
problems (in QCD subtraction of perturbation theory is a great one). Thre 
are two phases which are impossible in QCD: a case without the ground 
state (molecule force the system toward infinite Higgs VEV) and also a funny 
situation with chiral symmetry unbroken but confinement (hadrons exist but 
are degenerate in parity). It is amusing that chiral symmetry is restored at 
Nf = N c + 1, similar to what was found in QCD in the instanton model. 
Also note, that in this case the existence of the Coulomb phase is a proven 
fact. 



4 Inter-quark interaction and hadrons near 
and above T c 

It is well known that the nucleon effective mass in nuclear matter is reduced, 
and there are also indications that m p does the same R571 . Among approaches 
predicted "dropping masses" at high T the simplest (and the most radical) 
one is Brown-Rho scaling. According to it all hadronic masses get their scale 
from < qq >, and therefore vanish at T — > T c . This idea is supported by the 
QCD sum rules, provided the "vacuum dominance" approximation is used. 
However it is not so in a vacuum made of the instanton molecules, because 
they generate non-zero average values for some fermionic operators, even 
when T > T c and < qq >= 0. Those operators can be treated as new non- 
perturbative inter-quark interaction, which can be described by the following 
(Fierz symmetric, color-singlet only) Lagrangian [22 

I 



j mol sym 



G 



1 {^f-^T^f 



2N3 



(^r a 7 ^) 2 + $T a j^) 2 } + ^(^75^) 2 } + m) 



c c 



which has Numbu-Jona-Lasinio-type form and the coupling constant is pro- 
portional to density of molecules G = J n(pi, p2)dpidp 2 ^73-(27rp 1 ) 2 (27rp 2 ) 2 - 

I A 

Here, n(pi, p 2 ) is the tunneling probability for the IA pair and Tia is the corre- 
sponding overlap matrix element, r a is a four- vector with components (r, 1). 
The effective Lagrangian (fj) was determined by averaging over all possible 
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molecule orientations. Near the phase transition, molecules are polarized and 
all vector interactions are modified according to (ip^^Til)) 2 — > 4 (^70^) 2 . 

Numerical simulations and analytic studies have been used to 

calculate both spatial and temporal correlation function. The former expo- 
nentially decay with the so called "screening masses": their T-dependence 
for a number of hadronic channels is shown in Fig.(|4])(b) || show overall 
agreement with lattice ones. Especially important is strong attraction in 
scalar-pseudoscalar channels, shifting these masses down from their high-T 
asymptotic, M/ttT = 2. This attraction is also clearly seen in temporal cor- 
relation functions as well. For example, the pion correlators at several T are 
shown in Fig.(§)(a). They all have upward enhancement, which is absent in 
other channels (such as p). Fit of the correlator suggests that pion survives 
the phase transition and exist at T > T c as a (non-Goldstone) massive bound 
state. Most other hadrons "melt" into constituent quarks. Those have no ef- 
fective mass but effective energy (or "chiral mass") above T c . Unfortunately, 
accuracy of the calculations done so far does not allow to get information 
about the fate of vector and axial mesons, the major players in the next 
section. 

Unfortunately, the temporal correlators are very difficult to get on the 
lattice: there are only few points in time direction. The best data we can 
compare with definitely show attraction in the pion channel at T > T c , 
which is absent in vector ones. 



5 The U(1)a restoration? 

The fate of the U(1)a anomaly at finite temperature has received a lot of 
attention lately. A number of authors have emphasized that a (partial) U(1)a 
restoration may be the case, with rather dramatic observable consequences 
[II], |2"9"| . Indeed, recall that at T=0 the 7/ — -n mass difference is larger 
than all other meson mass splittings, and any tendency for it to shrink would 
therefore strongly affect the whole spectroscopy^. 

7 Those are directly analogous to the BCS theory of superconductivity, but with the 
instanton-induced interactions like (1). 

8 An interesting situation arises for Nf > 3 massless flavors. In this case, the anomaly 
does not affect the 77' correlation function above T c |3l], |2^, |3q ]. If chiral symmetry is 
unbroken, extra zero modes cannot be absorbed by the condensate, and the 't Hooft 
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It is well known that this splitting and the U(1)a anomaly are related with 
the topological charge, or well isolated instantons. It was argued in |J that 
those should be very rare in the ensemble at T > T c , because all instantons 
are "paired" into molecules. The density of isolated instantons above T c is 
small 0(m N f), but due to zero modes the quark propagators contribute a 
factor 1/m, and some results are finite in the chiral limit[]. So, a significant 
drop in the strength of the U(1)a anomaly around T c is expected, although 
to non-zero value. It was further suggested in that instead of dealing with 
7/ (which involve a complicated flavor mixing pattern), it is more convenient 
to look at a difference between ir and its U(1)a partner isovector scalar S. 

Such measurements have been done, both on the lattice and in the instan- 
ton model. All of them observe that the difference of it — 5 susceptibilities 
indeed drops dramatically at T c , indicating a move toward the U(1)a restora- 
tion. 

However, the calculations disagree about what exactly is its value at T > 
T c . The problem is the m — > limit, which is very difficult to do in practice. 
Then the conclusion is sensitive to extrapolations: MILC collaboration |32 



has fitted it quadratically x-k ~ Xs = C\ +m 2 C2 and gets a non-zero CI, but 
extensive results from the Columbia group (see their data in N.Christ's talk) 
have shown that the linear fit is much better, and it leads to Xtt — Xs ~" > 0- 

In the instanton calculations [[| the value is clearly non-zero (although 
it is not accurately determined yet). So, it seems that here we (for the 
first time) see a serious qualitative disagreement with lattice simulations at 
T > T c . This can be seen explicitly from the form of the Dirac spectrum. In 
the one we have calculated there are still quasi-zero modes due to isolated 
instantons, while in Columbia ensemble those seem to be completely absent. 
Why this happens remains unclear. 

The problem of rjr]' mixing at finite T was studied in more detail in 
[3C]. The most important conclusions is that there is a tendency towards 
ideal mixing, a separation of strange and non-strange component. However, 
contrary to others, it is shown that for T > T c the anomaly only can operate 
in the non-strange sector, with uu transition into dd being proportional to 
77V It is claimed that non-strange component may be rather heavy. 



vertex only contributes to 2iV/-point correlators. For Nf = 3 this means that the rj and 
7r are degenerate above T c , but the singlet and non-singlet A are split. 

9 Another way to explain it is to say, that for m=0 the individual instantons are absent 
in the vacuum, but can be created by the operators themselves. 
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In summary, a clear tendency toward diminishing of the role of U(l) 
anomaly at T T c is observed, but quantitative results for T > T c are still 
missing. If (as Columbia data suggests) it is vanishingly small, we would 
have 8 rather then 4 massless modes (for Nf — 2,m — 0), and understand 
why it moved away from 0(4) indices to the first-order ones. However, the 
instanton calculations show that it is not so small. 



6 The phase transition at finite baryon den- 
sity: the old and the new puzzles 

First attempts to introduce the non-zero chemical potential associated with 
baryon number have been done in quenched approximation a decade ago 
It was found that a very strange thing (referred as the "old puzzle") 



happens. At small \i as expected, nothing depends on it, till some threshold 
fi c is reached. But instead of finding a threshold to be around the constituent 
quark mass fi c ~ itin/3 0, they have found it at \i c ~ Those are not 

too different if quark bare mass is large, but if it goes to zero the latter value 
vanishes. That means that in the chiral limit some states with non-zero 
baryonic number exist (and can be excited) at arbitrarily small fi. There are 
no such states in the real world, so why lattice data have found them? 

For ten years it was not clear why, but the phenomenon was not a numeric 
artifact, and it has been since then seen in many other calculations. The 



resolution of the puzzle was recently made by Stephanov [35[ . He has pointed 
out, that it is "quenching" the QCD partition function which one should 
blame. Interestingly enough, it turns out that elimination of the quarks can 
be done in two very different ways. In both one writes in the bosonic partition 
function the factor [det(iD + ifijo + im] Nf coming from the integration over 
fermions, and then takes Nf — > 0. This factor is real at \i — and complex 
otherwise, so it can be written as a modulus and a phase. Stephanov has 
shown that the usual quenching corresponds to the limit, in which one keeps 
modulus and ignores the phaseQ. 

10 By tradition, lattice people just put /i in the Dirac operator without extra coefficients, 
which means that in their units quarks have baryonic charge f , not f/3. 

11 One can do this also without the limit, at finite Nf. such approach was suggested and 
studied separately by Gocksch p&. 
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Putting the absolute value of the determinant means that anti-quarks 
have a conjugate operator, with the opposite sign of the chemical potential 
fi — > —fi: so instead of having the baryonic charge -1 they also have it equal 
to 1, like quarks. Then mesons qTq may have the baryon charge 2, and 
the massless excitations excited at small /i are nothing else but such pions. 
Those exist in the Stephanov (unwanted) phase, in which the baryon number 
is spontaneously broken by the quark condensate[5 which also obtains the 
baryon charge. 

Another talk we had on the lattice data at finite baryon density have been 
given by M.Lombardo. This time it is based on rather complicated algorithm 
and strong coupling limit, but it is unquenched and includes the determinant 
as is, with its phase. 

A good news is the excitation to non-zero baryonic density happens 
around the expected threshold value. But the bad news is the excitation 
starts about 30% lower it and is more smooth compared to what one expects 
on the basis of the temperature value involved. So, we still see some un- 
usual baryonic states being excited, like some nuclei or nuclear matter with 
unexpectedly large binding energy. I do not know any explanation of this 
puzzle. 



7 "Dropping masses" in the experiment? 

Already at QM93, CERN dilepton experiments have provided first prelim- 
inary indications for a signal, significantly exceeding the theoretical expec- 
tations (e.g. @). Later CERES has found dramatic excess of dileptons at 
M e+e _ < m p . In spite of multiple attempts by theorists to explain it by 
"conventional sources" , it was found to be impossible []]. A list of "uncon- 
ventional" explanations (more or less in a chronological order) include: (i) 
dropping m p [ |3D| , pE0| ; (ii) pion occupation numbers at low momenta are very 
high [jH]| ; (iii) a very long-lived fireball [j38[| ; (iv) dropping [11 , |12[ ; (v) a 



12 In fact Stephanov has found boundaries of this phase and many interesting details 
about distributions of Dirac eigenvalues at non-zero [i, both in the wrong (no phase) and 
right (with the phase) ensembles, using a particular model of chiral restoration based on 
random matrix approach. Unfortunately, we have no time to discuss it here. 

13 In fact, I have never before seen that rather involved calculations by several groups 
agree so well, as far as the shape of spectra of M e+e _ is concerned. 
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modified pion dispersion curve (vi) dropping mA 1 

I will only discuss the p — A\ story, with brief remarks about others, (ii) 
The observed \ow-p t excess of pions (over thermal spectra) can indeed lead to 
a dramatic increase in low-mass dilepton yield. However, it is most probably 
due to resonance decays or effects of collective potentials. Both are late-stage 
phenomena, which hardly affect the early-stage dilepton production, (iii) 
Even if the long-lived fireball would appear at 200 GeV/A (which is hardly 
possible) it was found to produce about the same M e+e _ spectrum as the 
usual space-time scenarios, so it does not work, (iv) Above we have indeed 
suggested that m v /(T) should drop more than any other mass: however in 
order to account for CERES data one should increase yield of "escaping r/" 
by too huge a factor. (In connection to U(l) restoration issue, however, it 
would be extremely interesting to measure the rf yield, though.) 

The "dropping m p " idea is very well known, and it seems to be about 
the only one which can explain the dilepton data. It was studied in details 
by Li,Ko and Brown H(| in the cascade model based on Walecka-type model 
with attraction mediated by a scalar field. With much simpler hydro-based 
approach we have also verified, that one can get a very good description of 
CERES data by making the rho mass T-dependent, without any changes in 
standard thermal rate formulae or in space-time evolution. 

A point I want to make though is that "dropping m p (T)" are not con- 
sistent with chiral symmetry restoration, unless modifications of its chiral 
partner a\ follow. Relation between the two were made e.g. in the contents 



of Weinberg-type sum rules pH | . Both states should become identical at 
T c , so m p {T c ) = rriA^Tc). Important role of a\ for production of photons 



and dileptons was discussed in [0, f|I| [46]: we have used rather general 
expressions recently derived in |^7j Possible scenarios of how chiral restora- 
tion may proceed are shown in the m p — mA 1 plane in Fig.|5|(a). For example, 
path G corresponds to ref.|| , while C corresponds to Brown- Rho scaling. 
The corresponding dilepton spectra are shown in FigJ^(b): the variant D, 
with m p (T c ) = rriA^Tc) ~ (l/2)m p (0), does the best job, and Fig.|(c) we 
show contribution of separate stages in this scenario. The A\ contribution 
is mainly at small M e+e _, where CERES acceptance is low and background 



14 In order to explain why a\ is important, let us go "backward in time": it is the first 
hadronic resonance which may be excited in a collision of a photon, real or virtual, with 
a pion. 
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rather high. Still, with some effort one may be able to dig out the A\ contri- 
bution, using other information (especially the pt distribution). Note finally, 
that none of the scenarios mentioned above happen to violate current upper 
limit on the direct photons are set by WA80 pH . 



8 Flow, the "softness" of the EOS, and the 
fireball lifetime 

Our last topic deals with a more straightforward approach to observable 
signal of the QCD phase transition: instead of hunting for "dropping masses" 
we look instead at the effect of "dropping pressure". The ratio p/e, pressure 
normalized to the energy density, is shown in Fig.|6|(a) we show how this ratio 
depends on e. The existence of the "softest point" , the minimum of this ratio, 
is clearly seen. Our second point: in those (hydro relevant) plot the curves 
with and without baryon density look similar in many models: compare the 
two curves in Fig.|6]. 

Effect of the "softness" of EOS at T w T c on transverse radial flow was 
discussed for a long time |53], [54|[3 . More radical (and more controversial) 



idea was proposed in [49], where it was suggested that it can affect the 
longitudinal expansion and thus the global lifetime of the excited system in 
some energy window. In this workshop we had a detailed talk by D.Rischke, 
where he has shown that the effect of "softness" of EOS increases the lifetime 
at RHIC energies by about factor 2. 



It was shown in H9| that, due to softness of the EOS, there exists a 
window of collision energies in which the secondary acceleration of matter 
is completely impossible. Thus, if stopping does occur, a very long-lived 
fireball should be formed. 

To study this idea semi-quantitatively, a simple 1-fluid relativistic hydro- 
dynamic model was usecf^l Our lattice-based EOS shown (in unusual form) 
in Fig.^|(a) have smooth crossover in a narrow temperature interval AT ~ 5 

15 The issue was strongly oversimplified: in fact, transverse velocity is the integral over 
the whole history, and initial "softness" is only part of the story. Our hydro studies have 
shown, that the mean flow velocity is not sensitive to EOS at early stages, although there 
are modifications of the shape of its ditribution. 

16 The major uncertainties are in the initial conditions: we assumed that half of the 
total energy is stopped. 
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MeV at T c = 160 MeV the energy density jumps by about an order of mag- 
nitude. 

As the collision energy is scanned down, from 200 A GeV (SPS) to 10 
A GeV (AGS), we found three radically different scenarios: (I) At SPS one 
starts in the QGP phase, therefore longitudinal explosion quickly restore 
ultra-relativistic longitudinal motion, even if stopping a la Landau takes 
place, and expansion resembles the Bjorken picture, (see Fig.|7] (a)): (II) As 
the initial energy density passes the softest point, the QGP phase disappears 
and so does rapid longitudinal expansion. Instead we find a slow-burning 
fireball. For the heavy nuclei and initial conditions we discuss, the burning 
front moves mostly in the compressed longitudinal direction (see Fig.[F] (b)). 
The total lifetime of the fireball appears to be nearly as long as predicted 
for RHIC, but of course it is kind of upper bound on the effect. (Ill) At 
still lower energies burning process becomes more spherical, and ordinary 
expansion is developed. The main result, shown in Fig.|6|(b), is a significant 
peak in a lifetime. 

Remarkably enough, these completely different scenarios lead to not-so- 
different spectra! With little fitting of the initial conditions, they can be 
made consistent with available data. Thus, one-particle spectra alone is not 
enough, and details about the space-time picture has to come from analysis 
of flow and HBT. 

How one can test these unusual predictions and try to locate the "soft- 
est point" experimentally? First of all, it seems that it is above the AGS 
domainF], so one should rather scan downward at SPS. The following ideas 
are discussed in literature: (i) Look for maximal lifetime (or minimum of the 
HBT parameter A) ||49|| ; (ii) Look for the minimum of the "directed" flow 
in the collision plane j5(|; (iii) Look for the nearly isotropic distribution of 
dileptons, produced in the long-lived fireball pESfl . 
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Figure 3: Schematic phase diagram of zero T QCD (a) and supersymmetric 
QCD (b) as a function of the number of colors iV c and the number of flavors 
N f . 

Squares show where lattice calculations have found the infrared fixed point, 
dots are where the instanton ensemble is purely "molecular", with the un- 
broken chiral symmetry. 
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Figure 4: (a) Temporal correlation function for the pion (divided by the one 
corresponding to free massless quarks) versus distance (in fm) . Solid triangles 
and squares are for T = l.,1.13T c , open triangles, squares and hexagons 
for T.43, 0.6, 0.86T C , respectively, (b) T-dependence of the screening masses 
(given in units if 7rT) calculated from th|e IILM. Note that at chiral restoration 
the partners (a, 7r), (rho,ai) become identical. 
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Figure 5: (a) Possible scenarios of chiral restoration in m p — m ai plane, (b) 
the dilepton spectra corresponding to them (CERES data are shown with 
background from hadronic decays subtracted); (c) the most favorable case D 
is shown in details, with contribution from different stages. The dominant 
one is clearly the long-dashed one, or hadronic part of the mixed phase. 
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Figure 6: (a) The EOS in hydro-relevant coordinates for resonance gas 
and QGP. The solid line is for zero baryon number, the dashed line is for 
fib = 0.54GeV. The minimum is the "softest point" discussed in the text, 
(b) Solid line (and left scale) show the lifetime in the center, the dashed line 
(and the right scale) shows the space-time volume occupied by the mixed 
phase (T « T c ). 
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(a) 129 A GeV Au+Au (b) 30.9 A GeV Au+Au 




z (fm) z (fm) 

Figure 7: Space-time picture of AuAu collisions at SPS energies (a) and at 
the "softest point" (b). Solid lines correspond to fixed energy density while 
the short-dash lines are contours of fixed longitudinal velocity. Q,M,H stands 
for Quark, Mixed and Hadronic phases. Note qualitative difference between 
the two figures, as well as different time scales. 



24 



